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Abstract 

We study massless one-loop box integrals by treating the number of space-time 
dimensions D as a negative integer. We consider integrals with up to three kinematic 
scales (s, t and either zero or one ofF-shell legs) and with arbitrary powers of propa- 
gators. For box integrals with q kinematic scales (where g = 2 or 3) we immediately 
obtain a representation of the graph in terms of a finite sum of generalised hyperge- 
ometric functions with g — 1 variables, valid for general D. Because the power each 
propagator is raised to is treated as a parameter, these general expressions are use- 
ful in evaluating certain types of two-loop box integrals which are one-loop insertions 
to one- loop box graphs. We present general expressions for this particular class of 
two-loop graphs with one off-shell leg, and give explicit representations in terms of 
polylogarithms in the on-shell case. 
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1 Introduction 



Box integrals play an important role in the perturbative description of 2 ^ 2 scattering 
processes. Classic examples at one-loop include the scattering of light-by-light [|I| and the 
scattering of partons 0] . Recent improvements of experimental measurements demand even 
more precise theoretical predictions and there is significant interest in determining 2 — > 2 
cross sections at the two-loop order. To achieve this goal requires the evaluation of certain 
master two-loop graphs, such as the planar double-box graph [|, ^, or some one- loop box 
integrals with bubble insertions on one of the propagators. 

In 1987, Halliday and Ricotta p suggested a method of calculating loop integrals based 
on treating the number of space-time dimensions D as a negative integer. Because loop inte- 
grals are analytic in D (and also in the powers of the propagators), this is a valid procedure 
and, although the intermediate steps may be carried out in negative D (and in particular 
series expansions can be made), D remains a parameter of the calculation and can be taken 
to be positive after integration. The problem of loop integration is replaced by that of han- 
dling infinite series. This idea was neglected for some time until Suzuki and Schmidt started 
a more systematic application of the negative dimension method (NDIM) to a number of 
two-loop integrals |0, three- loop integrals one- loop tensor integrals |^ as well as the 
one-loop massive box integral for the scattering of light by light |Q. In this last paper, 
Suzuki and Schmidt discovered that as well as reproducing the known hypergeometric-series 
representations of Ref. 0], valid in particular kinematic regions, hypergeometric solutions 
valid in other kinematic domains are simultaneously obtained. Of course, all of these solu- 
tions are related by analytic continuation. However, it is easy to envisage integrals that yield 
hypergeometric functions where the analytic continuation formulae are not known a priori. 
In these cases, having series expansions directly available in all kinematic regions may be 
very useful. 

Recently, we have generalised this method to describe massive n-point one-loop graphs 



with general powers of the propagators and arbitrary dimension D |TT|]. For graphs with 
m mass scales, q external momentum scales and n legs, we have written down a template 
series solution with (m + g + n) summation indices, together with a linear system of (n + 1) 
constraints. The template solution is completely general, while the constraints can be read 
off the specific Feynman graph. By solving the system of constraints, we obtain many 
solutions with {m + q — 1) summation indices, each of which can be identified directly as 
a hypergeometric function in the appropriate convergence region. The full solution in a 
particular kinematic region is formed by adding the solutions that converge in that region. 
It turns out that by keeping the parameters general, it is easier to identify the regions of 
convergence of the hypergeometric series and, therefore, which hypergeometric functions 
to group together. This has the additional advantage of allowing a connection with the 



general tensor-reduction program based on integration by parts of Refs. [|12|, |13| where the 
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tensor integrals are linear combinations of scalar integrals with either higher dimension or 
propagators raised to higher powers. It is the goal of this paper to consider massless box 
integrals and to obtain expressions in terms of hypergeometric functions valid for general 
powers of the propagators and arbitrary dimension. 

Our paper is arranged as follows. In Sec. we show how NDIM can be applied to 
construct the template solutions for one-loop box integrals together with the linear system 
of constraints that relates the powers of the propagators in the loop integral to the summation 
variables. We give the expressions for the solutions in different kinematic regions for massless 
scalar box integrals with one off-shell leg and for the on-shell case in terms of hypergeometric 
functions of one or two variables. In both cases, D is arbitrary and the propagators are 
raised to arbitrary powers. As an application of the general formulae, in Sec. |^ we consider 
a particular class of two-loop box integrals which are one-loop box graphs with bubble 
insertions on one of the legs. We give general formulae for the scalar integrals with three 
powers of propagators set to unity and one propagator (corresponding to the place where the 
one-loop insertion is made) kept arbitrary. In this case, identities amongst hypergeometric 
functions can be used to simplify the general expressions. We show how to evaluate the 
hypergeometric functions in the on-shell case and, by making a series expansion in e = 
(4 — D)/2, give explicit expressions in terms of logarithms and polylogarithms for the relevant 
two-loop scalar integrals. Finally, our findings are summarised in Sec. ^. 



2 The general massless one-loop box integral 

The generic massless one-loop box integral in D-dimensional Minkowski space with loop 
momentum k is given by 

Jf (z.,, m) = J j,..^ _ j,.. ^ (2-1) 

where, as indicated in Fig. |l], the external momenta hi are all incoming so that Eti K = 
and the massless propagators have the form 

= k'^ + iO, 

Ai = (^k + Y,^k^ +iO ^7^1- (2.2) 

The external momentum scales are indicated with {Qf }. In our case they are the Mandelstam 
variables s = {ki + ^2)^, t = {k2 + k^)'^ and the external masses kf = Mf. In this paper 
we will focus on box integrals with at most one off-shell leg, so that we have kf = for 
i = 1, 2, 3, and fc| = M^. For standard integrals, the powers Ui to which each propagator is 
raised are usually unity. However, we wish to leave the powers as general as possible. Later 



3 



V, 



k4 



V4 



V2 



Figure 1: The one- loop box diagram. 

on we will use these general expressions to derive some results for two-loop box integrals 
with one-loop insertions on the propagators. 

We can rewrite Eq. (|2.1|) using Schwinger parameters Xj, so that 

I4 {j^u ^2, z^3, i^A] {QD) = J Vx j — ^ exp l^^XiAij , 



(2.3) 



where we have used the shorthand 



(4 1 poo 
r=i r(z^i) ^0 



(2.4) 



with 
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a = J2u,. (2.5) 

i=l 

Performing the Gaussian integral in a straightforward way we have the usual Minkowski- 
space result for massless integrals 

If (z/i, z/2, z/3, ^4; {Ql}) = I T^^:^ exp(Q/P), (2.6) 

where 

P = xi + X2 + X3 + X4, (2.7) 
while for box integrals with one off- shell leg (fc| = M^) 

Q = X1X3 s + 1 + X1X4 M^. (2.8) 

As usual, in the physical region t < and s > 0. 

To evaluate the integral further, we adopt the suggestion of Halliday and Ricotta 
and treat the number of dimensions D as a negative integer. This is valid because the loop 
integral is an analytic function of D. We follow the approach suggested by Suzuki and 
Schmidt 0-0 and detailed in |TT[ by viewing Eqs. ( |2.3| ) and ( |2.6| ) as existing in negative 



dimensions. We make a series expansion in Xj in both Eqs. ( |2.3| ) and ( |2.6| ). The role of 
having D < is that the power of V is now positive allowing a multinomial expansion. 
Following the notation of |TT|, we have 



I4 {1^1,1^2, 1^3, 1^4,3, t,M^ 



00 n 



vx j 



ni,...,n4=0 ' 



d'^k (xiAi)"^ {x^A^Y^ {x^A^Y 

1-^^0/2 



Pi,.^4=o qilq2W- Pil---P4l 

91, ■■■,93=0 

with the constraint 

D 

qi + q2 + q3+Pi+P2+P3 + P4 = (2.10) 

that ensures that the power of Q and V match up correctly. The integers Pi and are 
introduced in making the multinomial expansions of V and Q respectively. If more than one 
leg is off shell, then there will be additional terms in Q leading to more summation variables. 
Similarly, if we take the —>■ limit, this is the same as fixing = in Eq. (p.9|). 



The Xi are independent variables so that for the equality ( p.9|) to hold, the integrands 
themselves must be equal. Therefore, by selecting the coefficient of the powers of x^'^% where 
Ui = —Ui, on both sides of the equality we find 

..i,^or(i + gOr(i + g2)r(i + g3)lMr(i+por ' ^ ^ 

91,. ..,93=0 

subject to the system of constraints 

qi + q3+Pi = -1^1, 

q2+P2 = -^2, 

qi+P3 = -1^3, (2.12) 

q2 + q3 + PA = -^^4, 

qi + q2 + q3+Pi+P2+P3+P4. = -D/2. 

There are seven summation variables and five constraints so that two variables will be uncon- 
strained. The procedure for developing the solution for the loop integral further is detailed 



in Ref. |]TT|. Each of the fifteen solutions of the system is inserted into the template solu- 



tion ( 12.111 ). For example, solving with respect to the indices {qi,q2}, we find 

pi = Z/2 + Z/3 + ^^4 + g2 - ^/2, 
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P2 = -z/2 - q2, 

Pi = z/i + Z/2 + z/3 + gi - L'/2, 

?3 = -qi-q2 + D/2-Ui-U2-U3-U4, 



which is then apphed to ( 2.11| ). F functions that depend on the unconstrained variables gi 



and q2 are converted into Pochhammer symbols 

r(z + n) , , 

because they are the most suitable way to write generalized hypergeometric functions. De- 
noting this solution as /■t'?!''?^} and introducing the shorthand notation 

i^ij = + i^j, I'ijk = 1^1 + + Vk, (2-14) 

we have 

^ ^ r(i + f-a)r(i + .,23-f)r(i + .234-f) 

^ - f , gi + g2) (^^3, gi) (t^2, q2) {s/My it /My 

giilo (l + Z/123 - f , qi) (l + ^^234 - f , 92) 

The second line can be immediately identified as Appell's F2 function (see Eq. ( |A.4| )) while 
the apparently divergent F-function prefactor can be rewritten using the identity 

where the index i runs over all of the F functions in the numerator and denominator. This 
identity holds provided we treat D/2 as an integer, as we have already done in making the 
multinomial expansion. We see that 

^(A-«.) = ^, (2.17) 
i=i ^ 



which is generally true for all solutions and is independent of the Vi. Applying ( p.l6| ) to ( p. 151) 
we find that 



= (_i)f (^,2)#-^n-#)r(#-.,23)r(#-.234 



V {u^)V {u,)V {D - a) 
D D D s t \ , , 

X F2(a--,i/3,^^2,l + Z/l23-^,l + ^^234 --,^,^j. (2.18) 



Similarly, the other fourteen solutions are given by: 

r (.123 - f ) r (.,34 - f ) r (f - r (f - r (f - .,3) 
r(i/i)r(i/2)r(i/3)r(j/4)r(L>-(7) 

^/D D D ^ D ^ D s t \ 

X F2 (--.23,y-.i2,y -.34,1 + y -.123,1 + y -.234, 

fri.<?i} - / ^^f .f-...s4 rAf2A--i ^ (^^234 - f ) r (f - .123) r (f - .34) 
" ^ ^ ' ^ r(z.2)rMrp-a) 

^ / D D D s t \ 

X F2 (^.1, .3, — - .34, 1 + .123 ~ 2"' 2" ~ M2 ' M2 ) ' 

w} _ / ^-.,,3 / ..2\-'-^ ^ (^123 - f ) r (f - .12) r (f - 1/234) 
" ^"'^ ' ^^^^^ r(.i)r(.3)r(D-a) 

^ { D D D s t \ 

X F2 (^.4, y - .12, .2, 1 + 2" - .123, 1 + .234 ~ J^Jp^Jp) ^ 

-(-1)^^^ ' l^J r(.i)r(.3)r(.4)rp-a) 
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D D ^ D s M 



2' 



X -^2 .4 - .2, — - .12, .2, - i^34, 1 + — - .123, TT^-, 



_ / / .^2Af-'-^34 r (.134 - f ) r (f - .123) r (f - .34^ 



r(i/i)r(z/4)r(D-a) 

D D . D s M 



2' 



X i?2 .134 - 77' ^3, .2, TT - ^34, 1 + .123 - TT' 



r (.123 - f ) r (.2 - .4) r (f - .12) r (f - 1/23) 



a ( D D ^ ^ D s 

X 5i I .4, y - .23, y - .12, 1 - .2 + .4, 1 + y - .123, I , 

^ ^f-.£(^-f)r(f -^123)r(f -.134 



r(i/2)r(i/4)rp-a) 

a ( D ^ D ^ D s 

X -Si I (T- y,.i,.3,l + .134- y,l + .123 ~ Y' ~ J ' 
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^ ( D D D 

X F3 I Z/2, Z/3, — - 1/34, — - Z/12, 1 + Y - ^^14, — , — I , 



r{P2,g3} 



-1)- s~"''"-*r 



i/2 





X 52 




D 


D 






^1^134 - 


2"' 




r{pi>P3} 
-'4 


t{P2,P4} 
— -'4 


(s ^ t, 






' ^i) , 


t{P3,<12} 

-'4 


7-{P2,gi} 
— -'4 


(s ^ t, 




^ Z/4, t^2 ^ 




r{Pi.93} 
-'4 


r{P4,g3} 
— -'4 


(s ^ t, 




^ Z/4, 1^2 ^ 




7'{92,g3} 

-'4 


= it'''^ 


[s ^ t, 




^ 1^4, t^2 ^ 




7'{P3,g3} 
-'4 


T{p2,qs} 
— -'4 


{s ^ t, 


1^1 


^ Z^4, 1^2 ^ 


^ ^^3) • 



r (z/134 - f ) r (z/4 - 1^2) r (f - z/34) r (f - z/^ 
r(z/i)r(z/3)r(z/4)r(D-a) 



(2.19) 



The definitions of tlie functions F3, if2, "S*! and S2 are given in Sec. |A.1| together with a table 
of their regions of convergence. 

We divide the kinematic regions up as shown in Fig. ^ 
region I : > |s| + \t\, 

region 11(a) : \t\ > + \s\ and > 

region 11(b) : \t\ > + \s\ and |s| > (2.20) 

region Ill(a) : |s| > + \t\ and > |t|, 

region Ill(b) : |s| > + |t| and \t\ > 

and, applying the convergence criteria of Table |l] to each of the fifteen solutions, we find 
that they are distributed as follows: 



in region I 

Jf (z/i, z/2, 1/3, 1^4; s,t,M^) 
in region 11(a) 

Jf (z/i, z/2, 1/3, 1^4; s,t,M^) 
in region 11(b) 

I^il^l,l^2,J^3,^4;S,t,M'^) 

in region 111(a) 
in region 111(b) 

jf (z/i, z/2, z/3, 1^4; S,t,M^) 



r{<?i.i32} 



{Pl,P4} . 7-{P4,g2} 
"T J-A 



+ 1 



{pi,gi} 



r{P2,P4} _|_ j{p2,<ji} _|_ j{p4,g3} _|_ jUim} 



r{P2,P3} _|_ j{p2,g3} _|_ j{p4,g3} _|_ j{gi,g3} 



j{Pl.P3} _|_ j{P3,'J2} _|_ j{pim} _|_ j{g2,53} 



j{p2,P3} _|_ j{p3m} _|_ j{pi,g3} _|_ j{g2,93} 



(2.21) 
(2.22) 
(2.23) 
(2.24) 
(2.25) 




Figure 2: The kinematic regions for the one-loop box with one ofF-shell leg. The solid 
line shows the phase-space boundary \s\ + \t\ — M^, together with the reflections \s\ = 
\t\ + and |t| = |s| + M^. The reflections are relevant for the convergence properties of 
the hypergeometric functions which only involve the absolute values of ratios of the scales. 
The dashed lines show the boundaries Isl = and \t\ — M^. 
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Some solutions are convergent in more than one region. For example, I^^'^''^^' and I^'^^''^^' 
are convergent in both regions 11(a) and 11(b) while /j^^'^^J' is convergent in both 11(b) and 
Ill(b). We also see that in region 11(a), two of the solutions (^/j^^'^*-^ and I^^'^'''^^^ contain 
dangerous F functions when 1/2 = 2/4. These divergences indicate the region of a logarithmic 
analytic continuation and can be regulated by letting z/2 = z/4 + 6, canceling the divergence, 
and then setting 5—^0. Similarly, the two divergent contributions in region 11(b) 

and ll^^''^^^^ also cancel in this limit. 

We can perform several checks of these results. 



Analytic continuation 

The solutions in the different regions are related by analytic continuations of the hy- 



pergeometric functions (see for example the appendix of Ref. |]Tl 



The i/i = limit 

By pinching out one or more of the propagators (which corresponds to setting i/j = 0) 
we obtain results for triangle or bubble integrals (see Ref. [|lll)- For example, if we set 
= t's = 0, then any term containing l/F(z/2) or l/F(z/3) is eliminated. In fact, only 
five solutions survive, one in each group. In each case, the hypergeometric function 
collapses to unity and we obtain the expected result for the massless-bubble integral 
with off-shellness in each of the five kinematic regions thereby spanning the whole 
of phase space 

n 

(2.26) 



/2^(z/i,z/4;M2) = (M2)^~"""n>i,z/4) 



where we have defined, for future reference. 



n^(/i,/) 



F(/. + /.--f)F(f-/.)F(f-/.-) 



[2.27) 



The massless box: (i/i, 1^2, 1^3, 1^4, s, t, 0) 

The limit M ^ can be taken whenever the kinematic region allows it, that is to say, 
m regions 11(b) and Ill(b), where < |s|, < \t\. These two regions are related 
by the symmetry (s t, vi ^ 1/4, z/2 ^ 1/3), so we focus only on region 11(b). Only 
three of the solutions survive, and we have: 



if \s\ < \t\ 

/f (^l,^^2,^3,i^4;S,t,0) = 



{P2,'?3} 



M2=0 



{P4,g3} 



FU 



f ) F (f - Z/134) r (f - Z/123 

F (z/2) F (z/4) Y[p-o) 



D 



X 3-F2 ( 1^1, 1^3, 0- - 1 + ^^134 



D 



1 + 1^123 



D s 
2"'~t 



M2=0 
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+ -1 - s-'^'^H 



o o T - f ) r {U, - V,) r (f - V,,) V 



D 



T{u^)T{u2)T{u,)T{D-a) 



/ D D D s 

X 3F2 ( Z/4, — - Z/12, — - Z^23, 1 + Z/4 - Z/2, 1 + y - Z/123, 

, ,1. , r(z/i34-f)r(z/4-^^2)r(f -z/i4)r(f -^.34) 

^ ^ T {v{)T {v,)T {v,)T {D - a) 

D D D 



X 



iF2 I 1^2, — - «^14, — - «^34, 1 - Z/4 + Z/2, 1 + y - Z/134, -- ) . (2.28J 



Similarly, taking the same limit for solution (|2.25| ) in region Ill(b), we find the 

result valid when |s| > \t\, which is also obtained by applying the exchanges (s ^ t, 
Vi ^ 1/4, V2 ^ ^z) fo Eq. ( p.28| ). Note that we could have obtained the same result 



by returning to the template solution ( |2.1lD with the system of constraints ( |2.12| ) and, 



after setting ^3 = 0, solved the on-shell box directly. In this case, there are two external 
scales, s and t, so that there will be six summation variables (pi, . . . ,P4 and gi, ^2) and 
five constraints yielding six solutions, three of which converge when |s| < |t|, again 
yielding Eq. (1212^ ) . 

As before, there are apparent divergences in the F functions when U2 = ^^4 that must 
be regulated. This is straightforwardly achieved for particular values of the parameters 
by setting 1^2 = ^4 + ^ and making a Taylor expansion. 

The iyi = l limit: /f (1, 1, 1, 1; s, t, M^) 

If we set the propagator power equal to one, then all the groups ( |2.21| )-( |2.25| ) give the 
correct answer 



e2 r (1 - 2e) St 



-t)-2Fi( 1,-6,1- 



s 



+ {-s)-SF, { 1, -6, 1 - 6, ~) - {-MY'2F, (-e, 1, 1 



tj ' ' - - y ' ' St ^ 

(2.29) 

where u is defined hy s + t + u = and e = (4 — Z})/2. To obtain this result we have 
returned to the series representation of the hypergeometric function and manipulated 
the series by repeatedly summing with respect to one summation index to obtain an 
2-Fi function, applied identities to change the arguments of the 2-^1 and rewritten the 
2F1 as a series. Then we sum with respect to the other index, and repeat if necessary. 
Eventually all of the hypergeometric functions of two variables can be reduced to 2F1 
functions. 
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< — k4 




Figure 3: A one-loop insertion into a one-loop box diagram. 



3 Application to two-loop box graphs 



The general results for one-loop box graphs presented in the previous section may be 
applied to give analytic results for two-loop box integrals when there are one-loop insertions 
on one of the propagators. As is well known, the effect of such insertions is to modify the 
power to which that propagator is raised. For example, we consider the two-loop integral 
shown in Fig. 0, with off-shell legs 



1 



D/2 J ^/.5 ^A«6 : 



(3.1) 



where the Ai are independent of the second loop momentum / and are given by Eq. ( pl^ 
while 



B2 



r + iO 

{l + k + ki + k2 + k^f + iO. 



(3.2) 



The momentum flowing through the bubble is k + ki + k2 + k^ so that the result of the 
integration over / is (see Eq. ( p.26| )) 



«/2 fif fif 



it: 



h (Ai5,Ai6;^4) = n (/i5,/i6) Ai 



-Ms -Me 



(3.3) 



where 11^ is defined in Eq. ( |2.27| ). In this way, the overall power to which ^4 is raised to, in 
the two- loop diagram ( |3.1| ), is /i4 -|- /^s -|- /ie + /i? — Inserting Eq. ( p.3|) into ( p.l|) we find 



D 



Jf (/Xi, /i2, Ai3, {/i4, /i5, /^6, Z^?}; {QD) = n^(/^5, /^e) Ia [f^U /^2, /^S, /^4567 " y ; {<5i } j , (3-4) 

where /i4567 = /^4 + + /^e + yU7- Results for diagrams obtained by pinching out one of 
the propagators are obtained by setting the corresponding /Xj— i>0. For example, one of the 
boundary integrals of Ref. is obtained as the special case of (|3.4|) , with fi^ = fij = (see 
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Fig. ^ (a)). Similarly, the two-loop diagrams with one-loop insertions on the other three 
propagators are defined in an analogous way so that 

J? Ai2, yW3, /^4}, /^5, /^6, {QD) = n^(/^2, yUs) (^/^1234 - y , /^5, /^e, /i?; {Qj^}) , 



Jf (/ii, /i2, {/is, /i4, /i5, /ie}, ^^^] {Ql}) = n^(/i4, /i5) (^/il, /i2, /i3456 " ^ , {Ql}^ 



(3.5) 



where the notation is obvious. 



For box graphs with only one off-shell leg, the symmetry of the diagram reduces the 
number of distinct integrals to two: 

-^4^ ({Atl5/^2,/U3,/i4},/i5,/i6,/i7;S,^,^^) = (/U?, At6,At5, {/i4,/i3,/i2,/il};t,S,M2) , (3.6) 

•>^4^ (yUl, {/^2,/i3,/i4,/i5},/^6,At7; S,^,^^) = </f (yU?, /^G, {/"S, /i4, /i3, /i2}, /^i; S, M^) , (3.7) 

so that it is sufficient to consider diagrams with insertions on the third and fourth propaga- 
tors. In the on-shell limit (M^ ~^0), there is the further relation 

J? (/ii, Ai2, {/is, At4, /i5, /^e}, /^7; s, t) = jf (/i7, /ii, /i2, {/is, /^4, /i5, /^e}; ^, «) , (3.8) 
so that for the massless box we only need to consider insertions on a single propagator. 

3.1 One-loop insertions in the one-loop box with one ofF-shell leg 

In this section, we further specify the values of the propagators in the general forms for the 
one- loop box graphs of Eqs. (|2.21| )- (|2.25| ): we fix three of the propagator powers equal to one, 
while the fourth power is kept free. Because of the symmetry properties of the integral 
we need only to keep either 1/4 or 1/3 general. 



Jf (1, 1,1,^^4; M^) 

This limit is appropriate for two-loop diagrams such as that depicted in Fig. |^. We 
choose to work with the solutions in region I, given by Eq. ( p.21| ).p| Each of the four 
solutions is an Appell F2 function which can be represented as a double Eulerian 



^ Although we start from the solution for |s| + |t| < M^, the same expressions can be obtained starting 
from any of the other kinematic regions. 
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integral (see Eq. ( |A.12| )). However, for this choice of the parameters, the F2 functions 
simphfy (see Refs. [|11|, 0) and we find 

/f (1,1,1, z/4;s,t,M2) 

oxf-2-.4. 2 x-ir(3 + z/4-f)r(f -3)r(f -2-1/4 
^ ' \ ) V ; r (1/4) r (D - 3 - z/4) 

/ D D s s 

X Fi 1,2 + 1/4-77,1,4--,— 



X 



2'^'^ 2'M2'M2-t 

r(3-f)r(2 + .4-f)r(f-2)^r(f-i-.4 
r (z/4) r (D - 3 - z/4) 

ir(2 + z/4-f)r(f -3)r(f -1-Z/4 



+ (-1)- t--'^-""^ [M^ -t 



V (1/4) r (D - 3 - z/4) 



X 2F1 1,1,4 



2 ' M2 - t 



, n 3. 2 x-4r(3-f)r(f -2)r(f -2-Z/4 

X 2Fi(l,z/4,3 + z/4-^,^^^). (3.9) 

Note that the value of D plays no role in simplifying the hypergeometric functions 
and the result given here is for general D. The remaining hypergeometric functions 
can now be manipulated using standard identities and the one-dimensional integral 
representations given in Sec. [A.2| can be used for specific evaluations. At this stage, a 
series expansion in e = (4 — Z})/2 becomes necessary. 

jf (1,1,^^3, M^) 

Similarly, for the case where i/i = z/2 = z/4 = 1 and z/3 is kept general, we obtain 

r2 



/4^(l,l,z/3,l;s,t,M' 

^ (-i)# (Mf-'(M'-^)-(M--r' "^"'T(z;i;-V, 



f-2/ 2 o x-ir(3 + z/3-f)r(f -2-Z/3'' 



^ D St 

X 2^1 1,^^3,3 + z/3 



2' (M2-s)(M2-t) y 

+ (_l)f sf-2-3if-2-.3(j^^2_^_^^2-f ^^2_^\-3-l 
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X 



r(2 + .3-f)'r(f-2)r(f-i-.3 



r (z/3) r (/} - 3 - z/3 

X 2^^! ( 1, Z^3,3 + Z/3 



r (D - 3 - Z/3) 



2 ' M2 - t 

. ^ 2 . / 2\'ir(2 + z/3-f)r(f -2)r(f -2-Z/3 
^ ^ ^ ^ r(z.3)r(D-3-z/3) 

/ D D D s t \ , , 

X F2 1, 2,1, l-z/3,3 + i^3 ,— ,— . 3.10 

V 2 2 2 M2 M2y ^ ^ 

In this case, one F2 function does not reduce simply and it is necessary to resort to 
the two-dimensional integral representation of Eq. ( |A.12| ) for explicit evaluation. 



3.2 One-loop insertions in the massless one-loop box 

We can also attack the problem in the on-shell box. Here we set ui = 1/2 = I's = 1 and keep 
i>4 general, which is appropriate for diagrams such as those shown in Figs. ^ and ^. Insertions 
on the other legs are given by the symmetry properties of the integral (see Eqs. 



We therefore choose to work with the solution valid when |s| < |t| since that contains no 
r functions that are singular when ui = z/3. In every case, the 3F2 functions of Eq. ( |2.28|) 
reduce to 2-F1 functions and we find 

If (1,1,1, 3, t) = (-l)f 2j {2 tLAl 



X 2F1 (1,1,4 



r (z/4) r (D - 3 - 2/4) 

D _s_' 
2' t. 

2 



, ,0 n . r 2 + i/4-f r(z/4-l)r f-l-Z/4 



X 2^1(1, ^ - ^ 



s 



D 




2" 


- 1 




r 


h- 


1/4 


\ ^ 


D 
2 



r (1^4) r (£) - 3 - 1^4) 



^ ^ r(3-f)r(i-z/4)r(f -2 



2 



r (D - 3 - Z/4) 



. (3.11; 
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Figure 4: Two-loop box diagrams with pinched propagators. 

There is still an apparent divergence as z/4 — > 1 which can be easily removed by manipulating 
the hypergeometric functions using the well known analytic continuations to obtain 

n r(3 + i.4-f)r(f -2)r(2-f)r(f -2-Z/4 

/4^ (1,1,1, 1/4; S, t) = (-1)^ t^-2-.43-l^_ 27 V2 ) \ 2) \2 



r(i/4)r(D-3-z/4)r 3-f 



D D s + t 

2,- ^ 

2 2 



X 2Fi(l,--2,--l 



r(2 + z.4-f) 


V 




-0 


2 

r 




T{D - 


-3 


-z/4)r 


(3- 


f) 





+ (-l)T 

/ D s + t\ 
X (^1,1.4, (3.12) 

We have checked that the same result can be obtained by starting from the general solution 
valid for \t\ < \s\. In this case, we must regulate the singularity as z/3 — > z/i by setting 
^3 = ^1 + The singularity as 5— >0 is canceled by analytically continuing the and, 
after taking the 6—>-0 limit, we recover Eq. ( |3.12|) . 



3.2.1 Explicit evaluation of two-loop box integrals 

To give more explicit expressions requires a more precise knowledge of z/4. For the two- loop 
diagram shown in Fig. || the value of is given by /i4 + yUs + /xe + — -j = n — where n 
is an integer. The case n = 2 corresponds to the simplest case /X4 = /iy = and /is = /le = 1, 



shown in Fig. ^ (a). Substituting this value in the general expression ( |3.12| ) and restoring all 
overall factors, we find that the expression for this two-loop integral in D = 4 — 2e is given 
by (see Eq. 



J4^(l,l,l,{0,l,l,0};s,t) = (-t)-2^-^2Fifl,-e,l-e ' + ^ 



2se3^H' ' ' s 

K2 
2se 



+ (-^r'^ A2Fifl,e,l-e,^V (3.13) 
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where the constants Ki and K2 are given by 



r(l + 2e)r(l -e)^ 
(l-2e)r(l-3e) 

r(l + 2e)r(l -2e)r(l + e)r(l-e) 



(3.14) 



(3.15) 



(1 -2e)r(l -3e) 

Note that by starting off with the NDIM approach, we have not actually had to perform 
any integrations to reach this result or make any assumptions about the smallness of e. The 
hypergeometric functions have one-dimensional integral representations (see Eq. (|A.10|) ) and 
can be expanded around e = in terms of polylogarithms. The necessary integrals are easily 
done 

2F1 (1, -e, 1 - e, x) = 1 + e log(l - x) - U2 (x) - Lis (x) + O (e") (3. 16) 



iFi (1, e, 1 - e, x) = 1 - e log(l - x) + e' 



-2 Lis 



X 



X — 1 



U2 (x) 



2Li3 



X 



,x — 1 

where the polylogarithms are defined by 



Lis (x) + -log3(l-x) 



Li2 (x) 



z 



x<l 



and 



Lis (x) 



1 ^ log(-2) log(l — xz) 
z 



z 



For X > 1, the following analytic continuations should be used 



' 1 \ 1 TT 

Li2 (x ± zO) = — Li2 ( — I log^ X H ± ivr log x x > 1 

\x/ 2 3 



+ 0(e^), (3.17) 

(3.18) 
(3.19) 

(3.20) 



X < 1 



Li3(x±z0) = Li3(^) -^log3(x) + ylog(x)±z|log2(x) x>l. (3.21) 



Similarly, the integral with only one pinched propagator, fi^ = fi^ = fie = I and /xy = 0, 
shown in Fig. ^ (b) is given by 

s + r 



Fi l,-e,l-e. 



rf(l,l,l,{l,l,l,0};.,t) = i-t)-'';^^' 

-2. 3^2 ^ , , . S+t 



i-sy 



452^3 



Fi 1, 1 + e, 1 - e. 



(3.22) 



The series expansion for the first hypergeometric function is given by Eq. (p.l6|) while the 
second can be obtained from Eq. ( p.l7| ) by using Gauss's relation between contiguous hy- 
pergeometric functions 



(/?-a)(l-x) 2Fi(«, /?, 7, x)-(7-«) 2i^i(«-l, P, 7, x) + {-f-P) 2Fi(a, /5-1, 7, x) = 0, (3.23) 
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such that 

,F,(1.;3 + 1. 1 - X) = + ft 1 - <,.). (3.24) 

Finally, the scalar integral for the bubble insertion /i4 = /^s = /ie = A^7 = 1 shown in 
Fig. I is 

Jf(l,l,l,{l,l,l,l};.,t) = (-t)'^- 2il + e)!!^l^ ^ ~ ^' ^ 

+ (-^) 4 (1 + 26) .3^3 (1' 2 + e> 1 - ^> — ) • (3-25) 



Once again, the series expansion for the first hypergeometric function is given by Eq. ( p.l6|) 
while the second can be obtained from Eq. ( |3.17D by repeated use of Eq. ( ^.24[ ). 



4 Conclusions 



In this paper we have evaluated one-loop massless box integrals with arbitrary powers 
of the propagators and with up to one off-shell leg as combinations of hypergeometric func- 
tions. The method we have used, first suggested by Halliday and Ricotta, has its roots in 
the analytic properties of loop integrals and, in particular, the possibility of treating the 
space-time dimensions D as a negative integer in intermediate steps. In Ref. 0] we have 
developed a general strategy for evaluating one-loop integrals in NDIM and we have pointed 
out some subtleties that can occur in the application of the method. For the box integrals 
we have considered here, with q energy scales, we have expressed the final result as finite 
sums of hypergeometric functions with q — 1 variables, that converge in the appropriate 
kinematic regions. The general results for one off-shell leg in the kinematic regions specified 
by Eq. ( [^.20| ) are given in Eqs. ( 2.21 )-( 2.25 ). Similar expressions for the on-shell case are 



given in Eq. ( [^.2(j| ). We would like to point out that no integration was actually necessary 
in obtaining these results. 

All of these expressions are valid for arbitrary powers of the propagators and are therefore 
relevant to classes of multiloop graphs where there are (multiple) one-loop insertions on the 
propagators. We have studied how these expressions are relevant to this type of two-loop 
graph and, in particular, two-loop graphs with three powers of propagators set to unity 
and one propagator (corresponding to the place where the one-loop insertion is made) kept 
arbitrary. With this choice of parameters, identities amongst hypergeometric functions can 
be used to simplify the general expressions. Explicit results in terms of hypergeometric 
functions are given for the one off-shell case in Eqs. ( p.9| ) and ( |3.1CI| ). In the on-shell case, 
the two-loop scalar integrals reduce down to two Gaussian 2F1 functions. Up to this point 
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we have not actually had to perform any integrations explicitly or make a series expansion 
in e = (4 — D)/2. However, to write the hypergeometric integrals in terms of logarithms and 
polylogarithms it is necessary to use an integral representation and make the series expansion 
in e. For the 2-^1 functions, the integral representation is one- dimensional and the integrals 
are well known. Explicit results for the graphs of Figs. ^ and § are given in Eqs. (|3.13| ), 
(|;22D and (jpl). 



It is clear that NDIM is an extremely efficient way of solving one-loop integrals. Fur- 
thermore, as we have shown in this paper and as Suzuki and Schmidt P, 0] have previously 
shown, NDIM can help in evaluating multi-loop integrals where there are one-loop insertions 
on one or more of the propagators. Whether or not NDIM can provide some non-trivial 
results for multi-loop graphs is an interesting, but still open, question. 
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A Hypergeometric definitions and identities 

In Sec. |A.1| we give the definitions of the hypergeometric functions as a series together 



with their regions of convergence. Integral representations for the 2F1, Fi and F2 functions 
are given in Sec. |A.2| while identities for reducing the Fi and F2 functions to simpler functions 



are given in Sec. A. 3 



A.l Series representations 

The hypergeometric functions of one variable are sums of Pochhammer symbols over a single 
summation parameter m 

( ^ \ ^ (a,m)(Am) 

^"^^^"'^'"'"^ = £ (7,-) ^ ^^-'^ 

^ ( o o, , ^ ^ («,m)(/?,m)(/j^m) 

3F2 7,7, X) = E (^,^)(y,^) (A-2) 
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which are convergent when \x\ < 1. 



The hypergeometric functions of two variables can be written as sums over the integers 
m and n: Fj, i = 1, . . . ,4 are the Appell functions, H2 a Horn function and Si and S'2 
generahsed Kampe de Feriet functions: 



Fi {a,(3,(3\-f,x,y 
F2 {a,(3,(3\'^,i,x,y 
F3 (a, a , /3, 7, X, y 

H2 {a,p,-f,-f',5,x,y 
Si {a,a',(3,-f,5,x,y 
S2 (a, a', /5, 7, X, ?/ 



{a,m + n){f3,m){f3' ,n) x"^ 
(7,m + n) ml n\ 



E 



E 



{a,m + n){P,m){i3' ,n) 



~ (a, m) (a', n) m) (/?', n) x™ 

m! n! 



E ■ 

m,n=0 



E 



(7, m + n) 
a, m + n) (/?, m + n) x™ y^ 



°° (a, m — m)(7, ?T,)(7', n) X™ 

m! n! 



E 

m,n=0 



E 

?i,n=( 
00 

E 



(5, m) 

a, m + n) (a', m + n) (/?, m) x™ 



m,n=o (7,m + ra)(5,m) m! n! 

(a, m — n) (a', m — n) {P, n) {P', n) x™ 

m,n=0 



m! n\ 



(A.3) 
(A.4) 
(A.5) 
(A.6) 
(A.7) 
(A.8) 
(A.9) 



(7, m — n) 

These series converge according to the criteria collected in Table |l]. The domain of 



Function 


Convergence criteria 




X < 1, \y\ < 1 


F2,Si 


X + ?/ < 1 


Fa 


-y/|xj + ^J\y\ < 1 


H2, S2 


- X + l/\y\ > 1, X < 1, \y\ < 1 



Table 1: Convergence regions for some hypergeometric functions of two variables, 
convergence of the Appell and Horn functions are well known. That one for Si and S2 may 



be worked out using Horns general theory of convergence |15 
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A. 2 Integral representations 



Euler integral representations of 2F1, Fi and F2 are well known [|I^-[|I7| and we list the 
relevant formulae here. 



2F1 (a,/3,7,x) 



^^^^ X duu^~\i-uy-^-\i-uxy 



r(/?)r(7-/3) Jo 

Re{p) > 0, Re(7 - p) > 0. 



(A.IO) 



r(7) 



r (a) r (7 — a) Jo 

Re{a) > 0, Re(7 - a) > 0. 



(A.ll) 



F2 {a,p,P',-f,y,x,y) 



r (7) r (y ) 



r i/3) r i/3') r (7 - /3) r (7' - /?') 

1 rl 



X ( du I dvuf'-^v^'-^l - uy-"-\l - vy'-^'-\l - ux - vyy 

(A.12) 



/O JO 

Re(/3) > 0, Re{(3') > 0, Re(7 - /3) > 0, Re(7' - /?') > 0. 



A. 3 Identities amongst the hypergeometric functions 

The Fi and F2 functions have the following reduction formulae which leave a single remaining 
Euler integral at most [p!^ - [|T7 



y 



X 



F, (a, /3, {3 + /3\ X, y) = il- y)-%F, [a, + 
F2 (a, 13, /?', 7, a, x, = (1 - a - 7, x 

F2 («, /3, «, 7', X, = (1 - I3,a-I3, 7', Y 

F2 (a, /3, /?', /3, 7', x, y) = {l- xy%Fr [a, P', 7', ^) 

F2 («, /5, a, a, X, y) = (1 - - yy^'2Fi (^P, «, 

F2 (a, P, P\ «, P\ X, = (1 - - X - y)-^ 

F2{a,P,P\P,P\x,y) = {l-x-yy'^. 



y 



xy 



(l-xyi-y)^ 



(A.13) 
(A.14) 

(A.15) 
(A.16) 

(A.17) 

(A.18) 
(A.19) 



21 



References 

[1] R. Karplus and M. Neuman, Phys. Rev. 80 (1950) 380; Phys. Rev. 83 (1951) 776; 
B. De Tollis, Nuovo Cim. 32 (1964) 757; Nuovo Cim. 35 (1965) 1182; 
V. Constantini, B. De Tollis and G. Pistoni, Nuovo Cim. A2 (1971) 733. 

[2] R.K. Ellis and J.C. Sexton, Nucl. Phys. 269 (1986) 445. 

[3] V.A. Smirnov, "Analytical result for the dimensionally regularized massless on-shell 
double box" , hep-ph/9905323. 

[4] V.A. Smirnov and O.L. Veretin, "Analytical results for dimensionally regularized mass- 
less on-shell double boxes with arbitrary indices and numerators", hep-ph/9907385. 

[5] I.G. HaUiday and R.M. Ricotta, Phys. Lett. B193 (1987) 241; 

R.M. Ricotta, 'Negative dimensions in field theory', J.J. Giambiagi Festschrift, ed. 
H. Falomir, p. 350 (1990). 

[6] A.T. Suzuki and A.G.M. Schmidt, J. High Energy Phys. 9709 (1997) 002; Eur. Phys. 
J. C5 (1998) 175; "Solutions for a massless off-shell two loop three point vertex", 
hep-th/9712104; Phys. Rev. D58 (1998) 047701; J. Phys. A31 (1998) 8023. 

[7] A.T. Suzuki and A.G.M. Schmidt, "Negative dimensional approach for scalar two loop 
three point and three loop two point integrals", hep-th/9904195. 

[8] A.T. Suzuki and A.G.M. Schmidt, "Feynman integrals with tensorial structure in the 
negative dimensional integration scheme" , hep-th/9903076. 

[9] A.T. Suzuki and A.G.M. Schmidt, "Negative dimensional integration for massive four 
point functions. 1. The standard solutions", hep-th/9707187; "Negative dimensional 
integration for massive four point functions. II. New solutions", hep-th/9709167. 

[10] A.I. Davydychev, Proc. "Quarks-92", QCD161:A42:1992, World Scientific, Singapore, 
(1993), hep-ph/9307323. 

[11] C. Anastasiou, E. W. N. Glover and C. Oleari, "Scalar one-loop integrals using the 
negative-dimension approach" , hep-ph/9907494. 

[12] O.V. Tarasov, Phys. Rev. D54 (1996) 6479, Nucl. Phys. B480 (1996) 397. 

[13] O.V. Tarasov, Acta Phys. Polon. B29 (1998) 2655. 

J. Fleisher, F. Jegerlehner and O.V. Tarasov, "Algebraic reduction of one-loop Feynman 
graph amplitudes" , hep-ph/9907327. 

[14] A. Erdelyi, W. Magnus, F. Obcrhcttinger and F. Tricomi, "Higher Transcendental Func- 
tions", Vol. I, McGraw HiU (1953). 



22 



[15] H. Exton, "Multiple Hypergeometric Functions and Applications", Ellis Horwood 
(1976). 

[16] P. Appell and J. Kampe de Feriet, "Fonctions Hypergeometriques et Hyperspheriques, 
Polynomes D'Hermite", Gauthiers-Villars, Paris (1926). 

[17] W.N. Bailey, "Generalised Hypergeometric Series", Cambridge University Press (1966). 



23 



